Spin waves in the block checkerboard antiferromagnetic phase 
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Motivated by the discovery of new family 122 iron-based superconductors, we present the theo- 
retical results on the ground state phase diagram, spin wave and dynamic structure factor of the 
extended Ji — J2 Heisenberg model. In the reasonable physical parameter region of K2Fe4Se5 , we 
find the block checkerboard antiferromagnetic order phase is stable. There are two acoustic branches 
and six optical branches spin wave in the block checkerboard antiferromagnetic phase, which has 
analytic expression in the high symmetry points. To compare the further neutron scattering exper- 
iments, we discuss the saddlepoint structure in the magnetic excitation spectrum and calculate the 
predicted inelastic neutron scattering pattern based on linear spin wave theory. 
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I. INTRODUCTION 

Searching high-T^ superconductivity has been one of 
the central topics in condensed matter physical. Fol- 
lowing the discovery of the copper-based supercon- 
ductors two decades ago- , the second class of high- 
transition temperature superconductors has been re- 
ported in iron-based materials^""—. The iron pnictides 
contains four typical crystal structures, such as 1 Ill- 
type ReFeAsO (Re represents rare earth elements)^, 122 
type Ba(Ca)Fe2As2l'^, 111 type LiFeAa^ and 11 type 
FeSei^. The parent compounds with all structures ex- 
cept the 11 type have the stripe like antiferromagnetic 
state as the ground state. By substituting a few per- 
cent of O with F— or Ba with K^^, the compounds enter 
the superconducting (SC) phase from the SDW phase 
below Tc ■ In addition to the above four crystal struc- 
tures, a new family of iron-based superconducting mate- 
rials with 122 type crystal structure have recently been 
discovered with the transition temperature Tc as high 
as 33 K^^. However, these new "122" compounds differ 
from the iron superconductors with the other structures 
in many aspects^^. Firstly, the parent compound of the 
new "122" material is proposed to be Ko.8Pei.()Se2 with 
intrinsic root5byroot5Fe vacancy ordering determined by 
various experiments^^""— . Secondly, the ground state for 
KQ,sFei,QSe2 is Mott insulator with block antiferromag- 
netic order, which is observed by the neutron diffraction 
experiments^^""— . By the first principles calculation, Yan 
et. ali^ and Cao et. a& also found the block anti- 
ferromagnetic order is the most stable ground state for 
Ko,sFei,eAs2- 

To date, a number of researches have been carried out 
to study the nature of superconductivity and magnetic 
properties of these materials. The neutron scattering ex- 
periments by Bao et.al^^ have found the ground state of 
Ko,sFei,eAs2 to be block anti-ferromagnetic with mag- 
netic moment around 3.4/iB. It has been proposed by 
several authors that the magnetic and superconducting 
instabilities are strongly coupled together and the prop- 
erties of magnetic excitations, such as spin wave, play 



very crucial roles for the superconductivity in this family 
materials. Zhang et. aP^' even suggested that the su- 
perconducting pairing may be mediated by coherent spin 
wave excitations in these materials. 

In order to give the qualitative insight into the mag- 
netic excitation properties in this system, we studied the 
spin wave spectrum using the Heisenberg model on the 
2D square lattice with -^/S x -\/5 vacancy pattern. There 
are four independent parameters are used in the model, 
which correspond to the nearest neighbor and next near- 
est neighbor coupling between spins. We first obtain 
the ground state phase diagram as the function of those 
parameters, based on which we calculate the spin wave 
spectrum as well as the spin dynamic structure factor us- 
ing the Holstein-Primakov transformation. Our results 
demonstrate that the block checkerboard antiferromag- 
netic order is stable in a wide range of phase regime and 
there are two acoustic branches as well as six optical 
branches spin wave in this system, which can be mea- 
sured by the future neutron scattering experiments. 



II. MODEL AND METHOD 

The simplest model that captures the essential physics 
in Fe-vacancies ordered material K2Fe4^Se^ can be de- 
scribed by the extended Ji — J2 Heisenberg model on a 
quasi- two-dimensional lattice^ ■ ^ ^ , 
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Here, 6 = 1, 2, 3, 4 and 7 = 1,2, 3, 4; the first and second 
terms represent nearest-neighbor (n.n.) spin interactions 
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FIG. 1: (color online) Schematic diagram of the 2 dimensional 
crystal and magnetic structure (single layer for K2Fe4,Sez ) 
and the corresponding Brillouin Zone. On the left is the crys- 
tal structure and spin pattern of the pnictides K2Fe4,Se^ . 
The up-arrow(blue) and down-arrow(azure) atoms indicate 
the Fe-atoms with positive/negative magnetic moment, re- 
spectively. Here we show the considered model with nearest 
neighbor coupling Ji , J[ and next nearest neighbor coupling 
J2, J2. The coupling within each block is J(Ji, J2), and the 
coupling between blocks is J'(Ji, J2). On the right is the 
positive quadrant of the square-lattice Brillouin Zone show- 
ing wave-vectors F = (0,0), A = (^,0), B = (^, ^), 
C = (0, --^^). With the high symmetry line, the spin wave 
are contained along the direction T — A — B — T — C — B. The 
purple solid line marks the magnetic unit cell. The iron va- 
cancy site Fe is marked by the open square, and the occupied 
site Fe is marked by solid circle with the blue or reseda color 
indicating spin up or spin down. 

in the intra- and inter- block, respectively, as shown in 
FigUl The third and forth term are second-neighbor ( 
n.n.n.) spin interactions which are taken to be indepen- 
dent on the direction in the intra- and inter- block. Here, 
i is the block index, 7 denotes the nearest-neighbor block 
of i block. S'{S") represents the site-index which is n.n 
(n.n.n) site of site S. Ji and J[ ( J2 and ) indicate n.n. 
( n.n.n.) couplings of intra- and inter-block, respectively, 
which are illustrated in Fig. 1. Here, we define J2 is the 
energy unit. 

In order to understand this Ji — J2 Heisengerg Hamil- 
tonian, we depict the typical block spin ground state and 
the 2D Brillouin zone (BZ) in FiglTJ The q-vectors used 
for defining the high symmetry line of the spin wave is 
also depicted. Each magnetic unit cell contains eight Fe 
atoms and two Fe vacancy which are marked by the open 
squares, as shown in FiglTJ The block checkerboard an- 
tiferromagnetic order are recently observed by the neu- 
tron diffraction experiment in the K2FeiSe^ material. 
An convenient way to understand the antiferromagnetic 
structure of K2FeiSez is to consider the four paral- 
lel magnetic moments in one block as a supermoment; 
and the supermoments then form a simple chess-board 
nearest-neighbor antiferromagnetic order on a square lat- 
tice, as seen from FiglTJ. The distance between the 
nearest-neighbor iron is defined as a. Then, the crys- 



tal lattice constant for the magnetic unit cell is vTOa. 
FigHb is the 2D BZ for the magnetic unit cell. 

We use Holstein-PrimakofF (HP) bosons to investigate 
the spin wave of the block checkerboard antiferromag- 
netic ground state. As we know, linearized spin wave 
theory is a standard procedure to calculate the spin wave 
excitation spectrum and the zero-temperature dynamical 
structure factor— 1^. Firstly, we use HP bosons to replace 
the spin operators, as shown in Appendix A. 

Using Holstein-Primakoff transformations, we can ob- 
tain the HP boson Hamiltonian, 

^ = E V'l^fc^fe + E^-N- i?(fc^o) (2) 
fe 

Here, = (ali a+g a^^ &-fc-i fc-fc-2 b-fe-s b-fc-4); Eq 

is the classical ground state energy for block checker- 
board antiferromagnetic order, Eq = 8JiNS'^—4:J{NS'^+ 
4:J2NS^ - SJ^NS^. The Specific expression for Hk is 
shown in Appendix B. 

In the real space, we define the 'molecular orbital' , 
which is the combination of HP boson operators in one 
block. 

an = ^(oii— ai2+ai3 — ai4) 

ai2 = -^{aii+ai2—ai3 — ai4) 

an = -^{aii—ai2—ai3+an) 

an = -^(aii+ai2+ai3+ai4,) 

= ^(fc,-i-6i-2+6i-3-6i-4) 
Pi-2 = ^(fc,-l+b,-2-bi-3-6i-4) 
/3i_3 = ^(f),-i-6>-2-b,-3+6i-4) 

Pi-i = ^(6i-l+6,-2+6,-3 + f),~4) 

Here, the — represents the spin down block. The corre- 
sponding physical picture for each 'molecular orbital' is 
shown in the Fig. 2. 

In the 'molecular orbital' basis, the Hamiltonian be- 
comes, 

Hk = E^'fe^nr^^'^Vfe (3) 
fc 

Here, i^"^* = {al^ al^ a'l^ al^ fi.k-l l3-k-2 P-k-3 P-k-i) . 

The matrix elements for different 'molecular orbital' in 
the same block is zero, which is interesting; and we dis- 
cuss the Hamiltonian later. The Specific expression for 
^m-bitai shown in Appendix B. 

Because the boson Hamiltonian are big, one must use 
numerical diagonalization to solve eigenvalues for spin 
wave, which has a standard procedure^^"— , as shown in 
Appendix C. 

To further understand the physical properties in this 
spin system, wc obtain the analytical eigenvalues and 
eigenvectors in some special k points, such as fc = (0, 0) 
and k — ( ^ , ^ ) point. In the F point, the 'molecu- 

^ V lOa vlOa ' 

lar orbital' Hamiltonian becomes four 2x2 block matrix 
(/iJ!™+4 \Z**) about a™ and (m = 1,2,3,4), 

which indicates that the 'molecular orbital' are intrinsic 
vibration modes for F point. The spin waves at F point 
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FIG. 2: (color online) The schematic diagram for the 'molec- 
ular orbital' in the T point: (a) The deviation of spin in site 1 
and 3 has the same phase for the corresponding wave function; 
and the deviation of spin in site 2 and 4 has the same phase for 
the corresponding wave function. But the different between 
the phase of wave function for site 1 and 2 is 180 degrees. 
(b)The deviation of spin in site 1 and 2 has the same phase 
for the corresponding the wave function; and the deviation of 
spin in site 3 and 4 has the same phase for the corresponding 
wave function. But the different between the phase of wave 
function for site 1 and 3 is 180 degrees, (c) The deviation of 
spin in site 1 and 4 has the same phase for the corresponding 
wave function; and the deviation of spin in site 2 and 3 has 
the same phase for the corresponding wave function. But the 
different between the phase of wave function for site 1 and 2 
is 180 degrees, (d) The deviation of spin in site 1, 2, 3 and 4 
all has the same phase for the corresponding wave function. 



are collective excitations of the same 'molecular orbital' 
between different blocks. 

From the 2x2 block Hamiltonian, we obtain the eigen- 
values in the T point, 
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AFMl AFM2 AFM3 AFM4 

FIG. 3: (color online) The phase diagram for the Ji{J[) — 
J2(J2) model. The phases are defined in Ref. 17, among 
which the AFMl phase is the block checkerboard antiferro- 
magnetic phase observed in the neutron diffraction experi- 
ment. The blue/azure atoms indicate the Fe-atoms with pos- 
itive/negative magnetic moment, respectively. The magnetic 
configuration for AFM2, AFM3 and AFM4 is shown, (a) The 
phase diagram for the interaction parameter J2 = i, J2 ~ J2- 
(b) The phase diagram for the interaction parameter J2 = 
1, = 2.5 J2. 



gap in the T point. The final eigenvalue is also twofold 
degenerate. However, different from the above optical 
branches, it is a acoustic branch and always zero in F 
point, which is imposed by Goldstone's theorem. 

As shown above, we can also discuss the physical prop- 
erties in the k = {—^=-,—^=-) point. From the bo- 
son Hamiltonian, we can obtain the eigenvalues in the 

'= = (7k'7k)p°i"*' 



w'"'"' = -2SJi+S^(2Ji-j;-2.7^)2-,/{2 
w'"'"' = 2SJl+S^/{2Jl-J^-2J;^y-J^'^ 
w<''''''=cj<'"-''' = S-y/(2Ji-j;+2J2-24)2-4J^2_j/2 

There are also four Eigenvalues and each of them are 
twofold degenerate. Six eigenvalue are optical branches 
and two eigenvalue are acoustic branches. The third 
and forth eigenvalue are always degenerate in fc = 



-) point. 



^ VlOa ' VlOa ' 

Now, we have five different eigenvalues for spin wave 
in the special point, which can be used to fit the experi- 
mental data in order to get the n.n. ( n.n.n.) couplings 
of intra- and inter-block, Ji , J( ( J2 , J2 ) . Then using 
this interaction parameters, we can obtain the spin wave 
along all the BZ by numerical diagonalization method. 



There are four Eigenvalues. The first eigenvalue is 
twofold degenerate, and its eigenvector is a combina- 
tion by the first 'molecular orbital' in the spin up and 
spin down block site, as shown in Figl^ti- The second 
and third eigenvalues are degenerate and each of them 
is twofold degenerate. Their eigenvector are combination 
by the second and third 'molecular orbital' in the spin up 
and spin down block site, as shown in Figl^b and FigHJ;, 
respectively. They are all optical branches due to the 



III. RESULTS AND DISCUSSION 

In this section, we first present the phase diagram of 
the Ji — J2 model. Then we discuss the spin wave and 
spin dynamical factor in the block stipe antiferromag- 
netic phase. 

To investigate the phase diagram for the Ji — Ji Heisen- 
berg model, we use the stochastic Monte Carlo(MC) 
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method to investigate the system ground state. In the 
reasonable physical parameter region, the phase diagram 
for the Ji — J2 Heisenberg model is given in FiglSl which 
is plotted in the plane Ji/ J2 — J'l/ J2 at fixed value J2I J2- 
We obtain four stable phases in the case J2I J2 — 1- The 
first one is the block checkerboard antiferromagnetic or- 
der phase, denoted by AFMl in Fig.l, which has been 
observed by the experimenti^ii^. This phase is our inter- 
ested in the study of iron-based superconductors. Obvi- 
ously, when the coupling Ji is negative, the spin favors 
to form ferromagnetic configuration in the blocks, si- 
multaneously, when the coupling J( is positive, the spin 
favors to form anti-ferromagnetic configuration between 
the nearest-neighbor block. Of course, when J( is nega- 
tive, but small, the interaction J\ and J2 are dominant 
and the block checkerboard antiferromagnetic phase is 
also stable in this region. In the case J'21 J2 = 1, the 
block checkerboard antiferromagnetic phase stably exists 
when the following conditions are satisfied: Ji < and 
J{ > Ji . In the parameter region Ji > and Ji > J[ 
, the system favors to stay in the AFM2 phase, as seen 
in the Fig. [3l^a). In this phase, the antiferromagnetic 
order in the block arises from the antiferromagnetic cou- 
pling Ji. On the other hand, the spin system favors the 
AFM3 phase when J( < Ji < , which is mainly at- 
tributed to the dominant interaction J[ in this parameter 
regions. Similar to the phase AFM3, the system stay in 
the AFM4 phase in the region < Ji < J(, as seen in 
the Fig. [3ja). To concretely investigate the spin prop- 
erties in this system and compare the theory calculation 
with experimental results, in what follows, we change the 
parameter value J2/J2 = 1 to J2/J2 — 2.5 to investigate 
the phase diagram. In this parameter region, the antifer- 
romagnetic coupling J2 is dominant. Different from the 
first case, there are only two phase, AFMl and AFM2, 
in the phase diagram. The two phase is separated by 
the line Ji = 0.5J(. Below the line, the phase is AFMl 
phase, otherwise, the phase is AFM2 phase. It is interest- 
ing to ask in which region the realistic parameters of the 
iron pnictides fall. From the LDA calculations, Cao et al 
suggested that Ji = —29 mev, J[ — 10 mev, J2 — 39 mev 
and J2 — 95 mev^^. Such a set of parameters falls in the 
block checkerboard antiferromagnetic phase in FiglSl^a) 
and (b), implying that the K2FeiSe'^ should have the 
block checkerboard antiferromagnetic order. This fact 
tells us we only need to focus on parameter region in the 
AFMl phase. 

First of all, we use numerical diagonalization method 
to investigate the spin wave dispersion relations along 
the high symmetry direction in different situations. In 
the numerical calculation, it is convenient to set J2 = 1; 
and comparing with experiments, the actual energy scale 
of J2 for the specifical material can be deduced. Mo- 
tivated by the first principle reported parameters we 
first study the set of parameter: (i)Ji = —1, J[ ~ 0.2, 
J2 = 1 and J2 = 2.5. However, to study the influ- 
ence of interaction parameters on the spin-wave spec- 
tra, we also investigate four different sets of parame- 
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FIG. 4: (color online)Acoustic and optical spin-wave 
branches for linear spin- wave theory, T = 0, as a function 
of k along selected high symmetry directions in the magnetic 
unit cell Brillouin Zone. (i)Ji = —1, J[ = 0.2, J2 = 1 and 
J2 = 2.5, (ii)Ji = -1.5, J[ = 0.2, J2 = 1 and j!, = 2.5. , 
(iii) Ji = -1, J[ = 1.5, J2 = 1 and J2 = 2.5. , (iv) Ji = -1, 
J[ = 0.2, J2 = 0.2 and J2 = 2.5. , (v) Ji = -1, J[ = 0.2, 
J2 = 1 and J2 = 1.5. The x-axis correspond to the k point 
along the selected direction in Fig[T}3. y-axis correspond to 
the energy for spin wave; the energy unit is J2 = 1 except 
case (iv). 



ter: (ii)Ji = -1.5, J[ = 0.2, J2 = 1 and 2.5. 
, (iii) Ji = -1, J[ = 1.5, J2 = 1 and 4 = 2.5. , 
(iv) Ji = -1, J{ = 0.2, J2 = 0.2 and = 2.5. , (v) 
= -1, J[ = 0.2, J2 ^ 1 and = 1.5. In Fig. H we 
plot spin wave dispersions along high-symmetry direction 
r - ^ - B - F - C - B in the BZ for different interac- 
tion parameters. In all cases, there are one acoustic and 
three optical spin- wave branches (each of them is twofold 
degenerate). For the acoustic branches, the gap of spin 
wave in F point is always zero due to Goldstone's the- 
orem. In the first case, the vibration mode for acoustic 
branches is shown in Fig[2Jd). It is the collective excita- 
tion mode for the forth 'molecular orbital' in one block 
with the forth 'molecular orbital' in other blocks. And 
the relative phase in a 'molecular orbital' is only depend 
on the momentum k and independent on the interaction 
parameters. However, the relative phase between differ- 
ent block's 'molecular orbital' is dependent on the specific 
interaction parameters. The optical gap in the F point 
is dependent on the specific interaction parameters by 
contraries. As discussed in EqlH two of the three optical 
branches are degenerate at F point. Away from F point, 
the two degenerate optical branches split. For example, 
with the increasing of k along the TB direction, the opti- 
cal branch for the second 'molecular orbital' (Fig. [^Jj) is 
almost no change. In contrast, the optical branch for the 
vibration mode of third 'molecular orbital' (Figl^J:) has 
obvious dispersion, which can be clearly seen in FigH^. 
The reason that causes the splitting of spin wave can be 
attributed that the vibration mode of different 'molecular 
orbital' is not isotropic and dependent on the momentum 
k. Therefore, different vibration mode has different be- 
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FIG. 5: (color online) 3D spin wave dispersion for the pa- 
rameter set of Ji = —1, J[ — 0.2, J2 = 1 and J2 — 2.5 in the 
extended BZ for one block, as seen in Fig[l] a. The energy 
is in units of J2. (a) spin wave dispersion for first optical 
branch, (b) spin wave dispersion for second optical branch, 
(c) spin wave dispersion for third optical branch, (d) spin 
wave dispersion for Acoustic branch. The x-axis and y-axis 
correspond to 3kx — ky and kx + Sky direction respectively. 



havior in different momentum direction. Similar to the 
the above two optical branches, the third optical branch 
is related to the vibration mode of the first 'molecular 
orbital', which is a independent branch. 

With the change of interaction parameters, we can in- 
vestigate the influence of the interaction parameters on 
the spin wave, as shown in FigHl Firstly, in all cases, 
the acoustic branches is always zero in F point and there 
are always twofold degenerate in F and k = ( ^^^ , ^^^ ) 
point. In the following, we compare the spin wave dis- 
persion relation by the change of interaction parameters. 
Comparing with the Fig|4](i) and Fig|4](ii), we find that 
with the increasing of Ji, the spin wave dispersions be- 
come bigger in different k point, especially for the acous- 
tic branch in the k = ( ^ , ^ ) point. From the FigHl 

^ VlOa vlOa^ ^ — 

(i) and Fig|4] (iii), we observe that the amplitude of the 
second optical branch becomes large with the increasing 
of J(. It also change the energy for spin wave in dif- 
ferent k point, especially in the k = ( ^^^ , -^7^) point. 

With the increasing of J(, the amplitude of the second 
optical branch becomes more and more clear. With com- 
paring the FiglU^i) and FigU (iv), we observe that the 
interval for the first and second optical branch becomes 
smaller with the increasing of J2. The energy for spin 
wave is also changed in different k point, especially in 
the k — (-^7=^, ^/jQg^ ) point. Like all the above cases, the 
spin wave dispersion also has some change in the inten- 
sity with the increasing of Jj. However, the interval for 
the first and second optical branch becomes larger, which 
is one of the most obvious features for increasing J2, as 
seen in FigHJ^i) and Fig|l](v). 



FigE] shows the typical 3 dimension spin wave spec- 
trum in the extended BZ for one block of the first set of 
parameters in Fig.(|3^); and this plots provide a general 
qualitative overview. Regardless of the specific param- 
eter values, a common feature of the spin-wave disper- 
sion is that there are twofold degenerate in F and k = 
(--^=--, -^7=^) point and it has one zero branch in F point. 
One acoustic branch and three optical branch can also be 
seen, which is also a common feature in this system with 
the fixed saddlepoint's structure. Generally speaking, we 
can determine the interaction parameters by comparing 
the spin- wave gap at F and k = ( ^ , Z- ) point with 
the experimental data; then we plot the spin wave us- 
ing this set of parameters, which can be used to compare 
with the inelastic neutron scattering experiment. 

As we all know, the neutron scattering cross section is 
proportional to the dynamic structure factor S™{k,uj). 
To further guide neutron scattering experiment, we plot 
the expected neutron scattering intensity at different con- 
stant frequency cuts in k-space. The zero-temperature 
dynamic structure factor can be calculated by Holstein- 
Primakoff bosons. In the linear spin- wave approxima- 
tion, does not change the number of magnons, only 
contributing to the elastic part of the neutron scatter- 
ing intensity. However, S^{k) and S^{k) contribute to 
the inelastic neutron scattering intensity through single 
magnon excitations. The spin dynamical factor associ- 
ated with the spin-waves is given by the expression, 

/ Tn=±{l, 2,3,4} 

x5{Lj-ujf) (4) 

Here |0 > is the magnon vacuum state and |/ > de- 
notes the final state of the spin system with excitation 
energy. ^k„i is the m-th component of the eigenvector 
al\0>M. 

Fig. |6] shows our predictions for the intensity of the 
dynamical structure factor of the block checkerboard an- 
tiferromagnetic order(an untwinned case) as a function 
of frequency. At low energies (w — 0.5), four strongest 
quadrate diffraction peaks are visible, which come from 
the acoustic spin- wave branch. With the increasing of the 
cut frequency (w — 1), four strongest quadrate diffrac- 
tion peaks disperse outward toward, which is also in 
the acoustic spin-wave branch range. At intermediate 
frequency (w = 2), the dynamic structure factor becomes 
a chain ring shape diffraction peaks, which is a common 
results of acoustic and the third optical spin- wave branch. 
However, the intensity of diffraction peaks around the F 
point is very weak. At high frequency(u; — 3.5), the 
chain ring shape evolves to nine strong circle diffrac- 
tion peaks, which come from the two optical spin-wave 
branch, which is degenerate in the F point. But the inten- 
sity of diffraction peaks around the F point is also very 
weak. Different from the low energy situation, in the 
higher frequency(a; = 4.5), the nine strong circle diffrac- 
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block checkerboard antiferromagnetic state observed in 
the iron-based superconductors. There are two acoustic 
branches and six optical branches spin wave in the block 
checkerboard antiferromagnetic spin system, which are 
the combination by the 'molecular orbital' in the F point. 
Then, we discussed the saddlepoint structure in the mag- 
netic excitation spectrum, which can also be measured by 
neutron scattering experiments. The explicit analytical 
expressions for the spin- wave dispersion spectra at T and 
k — (-^=-, -^=-) have been given. Comparison with fu- 
ture inelastic neutron scattering studies, we can obtain 
the specific values of interaction parameters. We have 
also calculated the predicted inelastic neutron scattering 
pattern based on linear spin wave theory. In addition, we 
have also studied the specific influence of each interac- 
tion parameter on the spin wave and dynamic structure 
factor. Neutron scattering experiments about the spin 
wave and the behavior of spin wave at the proximity of 
a quantum critical point deserve further attention. 



Appendix A: Holstein-Primakoff Transformation 



FIG. 6: (Color online) Constant-energy cuts (untwinned) of 
the dynamic structure factor S^"{k,iu) for parameters: Ji — 
-1, J{ = 0.2, J2 = 1 and J2 = 2.5, (ii)Ji = -1.5. The x-axis 
and y-axis correspond to Sk^ — ky and kx + 3ky direction in 
range (— 7r,7r) respectively. The constant-energy cuts from top 
to down = 0.5 J2 , I.OJ2, 2.OJ2, 3.5J2, 5.0 J2 and 9.0 J2 for the 
all the energy range. 

tion peaks don't disperse outward toward, but inward to- 
ward; simultaneously, the middle circle diffraction peak 
around the F point becomes weaker and vanishes with 
the increasing of the cut energy. The diffraction peak 
at this energy cut come from the two optical spin-wave 
branches, which is degenerate in F the point. At high- 
est frequency(a; — 9.0), four strongest circle diffraction 
peaks are visible at the four corners of extended BZ for 
one block and four weaker circle diffraction peaks are vis- 
ible around the F point and the middle of the four bound- 
aries, which come from the first optical spin-wave branch. 
The results are similar for other sets of parameters that 
have same ground states, with the main difference be- 
ing that the energy cuts must be changed to obtain the 
similar spin dynamical factor patten. 



The Holstein-Primakoff transformation, 



5*1 = {y2S- alai)ai w V2Sai 
Sr = aU^2S-ala.)^V2Sal 
St = (S-a\a.^ (Al) 

for i = 1,2, 3, 4, which is occupied by spin up; S repre- 
sents the size of the magnetic moment for each iron, 
and 

^1, = bUi^2S-bl,b^.)^V2Sbl^ 
SZ, - {,j2S-blfi^,)b^,^V2Sb^, 
SU = [bl.,b^,-S) {A2) 

for i — 1,2, 3, 4, which is occupied by spin down. 
The fourier transform for bosonic operators is. 



IV. DISCUSSIONS AND CONCLUSIONS 



In conclusion, starting with the Ji — J2 Heisenberg = N^YikC *'^'^^a|, 

Hamiltonian model, we have obtained the magnetic ^ _ N^Y, e^^'^^b 

ground state phase diagram by MQ approach and found ^ 1 _ ^ ^ 

that the block checkerboard antiferromagnetic order is "l ~ ^fe^ "fc 

(A3) 

stable at reasonable physical parameter region. In this 

paper, we have used spin wave theory to investigate Here, we define N is the number of the magnetic unit 

the spin wave and dynamic structure factor for the cell. L is the site index and k is momentum index. 
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Appendix B: Spin Wave Hamiltonian 

Wc write the detailed expression in Eq. 2. Replacing 
the spin operators by HP bosons, we can get the HP 
boson Hamiltonian "H^, 



The matrix elements for different 'molecular orbital' in 
the same block is zero. 



Appendix C: Numerical Diagonalization Method 



El Fk 


Gk 


Sk 





{Bk)* 


{Ck)* 


{Ak) 


{FkT El 


Sk 


Tk 


Bk 





{Ak)* 


{Dk) 


{GkT (SkT 
{SkT {TkT 


El 


{Fk)* 


Ck 
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Dk 
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Fk 


Gk 
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{Dk)* 


(Fk)* 


El 


Sk 


Tk 


Ck Ak 





Bk 


{Gk)* 


{Sk)* 


El 


{Fk) 


Ak Dk 


{BkT 





{Sk)' 


{Tk)* 


Fk 


El 



and 



Au 
Bk 
Ck 
Dk 
Fk 
Gk 
Sk 
Tk 



Jl + lj(-lJ2+272 

5Jjg-rf-(0.5il-0.555) 

-VOAat + VTlia~f 



In the 'molecular orbital' basis, the Hamiltonian ■^o'"'"*'*' 
becomes, 



and 



El 





Q Q Ao Bo 


a 







El 


{Bo)* Fo 




Go 








El {Co)* -Do 




So 








Et Do {GoT {SoT 


To 


At 


Bo 


Co Do El 








{Bo)* 


Fo 


-Do Go El 








{Co)* 


-Do 


Po So 


El 





Do 


{Go)* 


{So)* To 





E'o 




K = 


-16Ji+4j(+8.72 








El = 


-8Ji+4jJ-8J2+8J^ 








K = 


-8Ji+4J(-8J2+8J^ 








Et = 


4J(+8J2 








Ac = 


[2J[-4j!^].[cos(ao)+cos(bo)] 








Bo = 


— 2J-^^[s^n{ao)-\-s^n{bo)]^ 










+4J2Sm(6o)i 








Co = 


— 2J[ [sin{ao)~sin(bo)]i 










+4J2Sm(ao)i 








Do = 


2.J[ [cos (a o)^ cos {b o)] 








Fo = 


-2j[[cos{ao)+cos(bo)] 










— 4J2 [cos (a 0) — cos {b 0)] 








Go = 


2J[ [sin{ao) — sin{bo)]i 










+4J2sin(ao)'i 








Po = 


4J2 [cos(ao) — cos{bo)]] 










— 2.J[ [cos{ao)-\-cos{bo) 








So = 


AJ2Sin{ba)i 










+2 J( [sin{a ) + svii{b a)]i 








To = 


[2j[+ij'2].[cos(ao)+cos{bo)] 










1.5 k •xt.+Q.h k -yt 








■ = 


0.5 k -xt — 1.5 fc -y^ 







We can use the numerical method to diagonalize the 
boson pairing Hamiltonian. The matrix form for boson 
Hamiltonian is. 



Here, = {a^,b) = (^a\,al,al,a\,b-i,b-2,b-3,b-i^, 
A and B is a 4 x 4 matrix; a and b are boson operators. 
The operators satisfy the commutation relation. 



'A B' 




a 


B A 




_6t_ 



and 



/ 
-/ 



7 is a 4 X 4 identity matrix. 

The diagonalization problem amounts to finding a 
transformation T, which lets T'^T-LT become a diagonal- 
ization matrix Vt. 



a 




a 


_b\ 


= f 





Tip 



Here, we require that a and /3 are also a 

set of boson operators and ip"^ = {oi\l3) = 

{a\,a\, 4, a\, /3_2, /3-3, P-i) ■ 
Then, 



^Fi,i'I-i',j>Ty. 



we get 



f /_f t = L 



(Cl) 



Since P = L we know / = {Lf){Lf^). Due 
to {I-T) and {I-T"^) are each inverses, thus I = 
(/_f t)(J_f ) and ft/.f = L. 

Then we want the final form of H is diagonalization. 



8 



Then 



[a\b\ 



'A B' 




a 


B A 

















'A 
B 


B' 

A 


f 




Vt 


' A 

B 


B 
A 


[a\p] {h} 


a 



T 



a 



4 



(C2) 



Here,f2 = diag (wi, • • • , ^4, u-i, • • • , uj-4) represents a 
diagonalization matrix and the matrix elements for diag- 
onaUzation matrix is (wi, • • • ^lo 4^,10^1, ■ ■ ■ ,^-4) . 

In other words, we want the expression f^Uf = n and 
matrix f2 is a diagonahzation matrix. We must solve the 
matrix T. 

Combining Eq ICll and Eq |C21 we get, 



TI 



= n 



HT 



I-HT 
{i-H)f 

ii-H)f 



= T 



Here A — 



T X 



(C3) 



= diag (wi 

is a diagonalization matrix. In other words, if we want 
to get T^HT — Cl, we can solve the general Hamiltonian 



{I-'H)T = T 



J.L. van Hemmen's^i strategy is that the canoni- 
cal transformation T is fully determined by its n(=8) 
columns {xi, ... , xi, ... , X-4}. We, therefore, 

reduce Eq |C3l to an eigenvalue problem for these n(=8) 
vectors Xi, 1 < i < 8. 



(i^n) X = Ax 



Where x = Xj, 1 ^ i ^ 8, and A £ 
{a;i,W2,W3,W4,-a;_i,-a;_2,-a;_3,-w_4}. So, to diag- 
onalize the Hamiltonian is equivalence to solve a general 
eigenvalue problem. 

/ 



If we define a matrix ly ~ 



I 



identity matrix, it is easy to proof , 



and / is a 4 X 4 



I^l_4 



I—lk. — ly ^I—'fily 

ere is a eigenvalue Xi 
for the general Hamiltonian, 



We assume there is a eigenvalue A^ and eigenvector 



Then the x'l = ly ■ Xi = 



Ml-4 



is also eigenvector of 



I-H and the corresponding eigenvalue is A' 



-A: 



-iy^Ln 



I-U 



I-Hx 
iyXr 

lyXt 

lyXt 



\Xi 
KXi 

— lyXiXi 
-KlyXi 

-A.X^ 

Kx', 



Therefore, the eigenvalue are in pairs. For the sake of 
convenience, we arrange the order of eigenvalues by the 
relative size of the value =| /ii P — | Vi P, (Hi > 

> Hg); for the same H^, we arrange the order of 

eigenvalues by its relative size. If the eigenvalue Xi for the 
corresponding eigenvector Hi > H2 > H3 > H4 is lower 
than zero, it represents the ground state is not stable. For 



the first four eigenvector | Hi 



and for the last four eigenvector | /ii P 
(1=5,6,7,8). 



Vi \'= 1(1=1,2,3,4) 



V,. V= -1 
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